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A density functional theory study of substitutional carbon impurities in ZnO has been performed, 
using both the generalized gradient approximation (GGA) and a hybrid functional (HSE06) as 
exchange-correlation functional. It is found that the non-spinpolarized Czn impurity is under almost 
all conditions thermodynamically more stable than the Co impurity which has a magnetic moment 
of 2/iB, with the exception of very O-poor and C-rich conditions. This explains the experimental 
difficulties in sample preparation in order to realize d°-ferromagnetism in C-doped ZnO. From GGA 
calculations with large 96-atom supercells, we conclude that two Co-Co impurities in ZnO interact 
ferromagnetically, but the interaction is found to be short-ranged and anisotropic, much stronger 
within the hexagonal afe-plane of wurtzite ZnO than along the c-axis. This layered ferromagnetism 
is attributed to the anisotropy of the dispersion of carbon impurity bands near the Fermi level for 
Co impurities in ZnO. From the calculated results, we derive that a Co concentration between 2% 
and 6% should be optimal to achieve d°-ferromagnetism in C-doped ZnO. 

PACS numbers: 71. 15. Mb, 75.50.Pp, 61.72.U- 



I. INTRODUCTION 

Ferromagnetism is a quantum phenomenon 
which arises from the long range ordering of in- 
teracting magnetic moments in a solid. The rich- 
ness of magnetic properties in solids is due to hy- 
bridization of d- and /-orbitals which is decisive 
for magnetic moment formation and the way in- 
dividual moments interact with each other in the 
solid. This is one of the reasons why magnetic 
semiconductors or diluted magnetic semiconduc- 
tors (DMS) are traditionally conceived as semi- 
conducting materials with transition metal ions 
as dopant. Numerous studies of DMS, both ex- 
perimental and theoretical (see, e.g. Refs. [thill) 
have come up with varying conclusions. In ZnO, 
substituting the cation by an isovalent transition 
metal, e.g. cobalt, introduces magnetic moments, 
but in itself does not lead to ferromagnetic cou- 
pling between these moments @. Hence, complex 
strategies, for instance co-doping, and elaborate 
growth techniques seem to be necessary to fabri- 
cate a DMS based on ZnO as host material. 

Recently there are reports on magnetism in 
semiconducting materials without any d- or /- 
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elements p-Q. This is commonly known as d°- 
magnctism and discussed for different classes of 
materials, starting from pure carbon based ma- 
terials like graphite, graphene, and non-metallic 
nanoparticles [10|, LLl| ■ The common feature in this 
class of materials is that magnetism involves the 
2p orbitals. The origin of magnetic moment is at- 
tributed to incomplete cancellation of majority and 
minority spin electronic contribution due to the 
Hund's type exchange. Magnetic moments have 
been reported for almost all ionic semiconductors 
in the presence of intrinsic defects, such as cation 
vacancies, or in the presence of impurities related 
to anion substitution. In the following we discuss 
reports of magnetism in semiconductors with im- 
purities at the anion site. 

The formation of a localized magnetic moment, 
which is a pre-requisite for d°-magnctism, depends 
on the relative strength of electronegativity be- 
tween the dopant element and the anion of the host 
semiconductor. If the bond between the dopant 
and the cation is weaker than the native bond of 
the semiconductor, this leads to localized atomic- 
like 2p orbitals of the dopant and a stable spin- 
polarized state. On the other hand, if the strength 
of the bond is stronger than the native bond of 
the semiconductor, this will cause derealization of 
the dopant 2p orbitals due to strong hybridization 
with the cations. Consequently, there is reduced 
or vanishing spin-polarization in the system 
Hence, substituting the anion site by an element 
with smaller electronegativity can introduce spin- 
polarized states in the gap which can then mcdi- 



ate the magnetic interaction through the double- 
exchange mechanism as discussed by Mavropoulos 
et al. [13j | . 

By means of DFT calculations, Peng et al. flij 
have shown that many elements in the second row 
of the periodic table (X = B, N, C) when substi- 
tuted at the anion site can develop a spin polarized 
solution in A1N (Xn) and ZnO (Xo). However, 
the nature of the magnetic correlation in case of 
B and N is still under debate [l4l - [l7j] . Specif- 
ically for ZnO, Adcagbo et al. found that both 
No and Nz n develop spin-polarization with almost 
1 ^b/N in ZnO JT^] . Lyons et al. have stud- 
ied No in ZnO by hybrid-functional calculations 
and observed that No defects create a deep impu- 
rity state [l9|. These states, if partially filled, can 
be spin-polarized. Also for Co, DFT calculations 
predicted a magnetic moment for the neutral or 
singly charged impurity [2(j • There are experimen- 
tal reports of room temperature ferromagnetism in 
carbon-doped ZnO films grown by pulsed-laser de- 
position [2l|, [22j ■ Measurements of the anomalous 
Hall effect have shown that the ferromagnetism is 
mediated by the charge carriers in ZnO. Remark- 
ably, results have been reported both for samples 
showing n-type conductivity 21[ or p-type conduc- 
tivity |23[. Moreover, Yi et al. reported on room- 
temperature ferromagnetism in carbon doped ZnO 
films when co-doped with nitrogen [24[ . 

In this work, we select ZnO as a host material 
for the study of d° magnetism. Impurities and 
defects in ZnO have been an active topic of re- 
search [25|, [26|. Here we focus on the magnetic 
properties of ZnO with carbon as substitutional 
impurity (on the zinc site denoted by Cz n , and 
on the oxygen site denoted by Co) studied with 
the help of density functional theory (DFT) cal- 
culations. The paper is arranged as follows: In 
the next Section we present technical details of 
our DFT calculations. The Results and Discus- 
sion Section discusses the stability of carbon im- 
purities in ZnO by calculating their formation en- 
ergies. The distance-dependent exchange interac- 
tions of Co in ZnO are obtained from the total 
energy differences of ferromagnetic and antifcrro- 
magnetic states. Finally, we conclude on the pos- 
sible connection between C impurities and the fer- 
romagnetism in ZnO found in experiments. 



II. 



COMPUTATIONAL DETAILS 



Electronic structure calculations are performed 
using the plane-wave pseudopotential DFT 
method as implemented in the Vienna ab-initio 
simulation package (VASP) [13, H^. For studying 
a substitutional impurity in ZnO, the system 
is modeled as a supcrcell consisting of periodic 




Figure 1: (Color online) Supercells and impurity con- 
figurations used in the present study. The small (blue) 
spheres represent Zn, and among the large spheres, the 
darker ones represent carbon on the O sites while light 
spheres represent oxygen atoms in the ZnO lattice. 
The supercells S122 (a), S221 (b) and S222 (c) contain 
a single impurity and are used to model impurity con- 
centrations of 12.5% for S122 and S221 and 6.25% for 
S222. The Co-Co interaction is studied with two Co 
in S226 (d) and S622 (e). One of them fixed at the po- 
sition CO and AO in S226 and S622, respectively, and 
another placed on one of the C1-C13/A1-A13 sites, 
which corresponds to a concentration of 4.17%. 



repetitions of the primitive wurtzitc cell. Since 
we use different sizes of supercells, we stick to a 
particular nomenclature in our discussions. We 
refer to the size of the supcrcell with periodic 
extention n x xn y xn z as "Sn.j;n y n 2 " . The elec- 
tronic structure of homogeneously distributed 
impurites has been studied with the relatively 
small supercell S222, see Fig. 1(c), while to study 
impurity pairs in ZnO we employ two kinds of 
large supercells, namely S622 and S226 (Fig. 1(d) 
and (e), respectively. This is necessary as the 
wurtzite crystal structure is non-centrosymmetric 
and the bond lengths are different along the 
hexagonal plane and along the c-direction. Es- 
sentially, S622 intends to scan the interactions 
along the hexagonal plane of the lattice and the 
S226 scans the interaction along the c-axis of the 
lattice. 

The calculations have been performed us- 
ing a plane-wave energy cutoff of 400 eV. We 
have used the generalized gradient approxima- 
tion (GGA) with the PW91 parameterization [29| 
for the exchange-correlation potential and the 
ion-electron interactions are treated with the 
projector-augmented wave (PAW) method f3Cj . 
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The lattice constants of the primitive lattice of 
ZnO is taken to be a = 3.29 A, c/a = 1.606 and 
the relative shift of Zn and O planes, u = 0.380, 
which were found to be the minimum e nerg y lat- 
tice parameters in the GGA calculations [32| • The 
doped supercell configuration is relaxed to avoid 
any strain until the structural energy is converged 
to 10 -6 cV. At this level of tolerance, the forces in 
the system are found to be below 3x 10 -4 cV A -1 . 
The relaxation is performed with a T-centered 
Monkhorst-Pack /c-points grid of 3x5x5 and 5x 
5x3 for supercells S622 and S226, respectively. 

Part of these results have been cross-checked 
with the full-potential local orbital (FPLO) 
minimum-basis code, which uses localized basis 
sets [H|. DFT within GGA gives reliable struc- 
tural parameters for ZnO but has limitation in pre- 
dicting the right optical band gap, which is with 
a value of 0.69 eV grossly underestimated in com- 
parison to the experimental value of 3.4 eV [32| . 
In addition, the localization of Zn d-bands is not 
properly described by GGA because the hybridiza- 
tion of Zn d and O p bands is overestimated. Thus, 
the Zn d bands are located almost 3 eV too close 
to the Fermi level as compared to experimental 
observation (33[. One way to overcome the prob- 
lem is through the introduction of hybrid func- 
tionals where the electronic exchange potential is 
mixed with some percentage of Hartrcc-Fock ex- 
change, while the electronic correlation potential 
is entirely taken from the GGA. We have em- 
ployed the screened hybrid-functionals to ascer- 
tain the electronic structure of a single impurity. 
The hybrid-functionals are of the type suggested 
by Hcyd, Scuseria and Ernzerhof [34| , with the the 
exact (Hartree-Fock) exchange and the GGA ex- 
change in the ratio 1:3 [35[. We have used the 
HSE06 functional for our studies (36j . 



III. RESULTS AND DISCUSSION 

The stability of an impurity configuration can 
be judged from its formation energy. Since substi- 
tution changes the number of atomic species, the 
formation energy must be calculated with reference 
to thermodynamic reservoirs characterized by their 
chemical potentials. In practice, the values of these 
chemical potentials reflect the experimental condi- 
tions during sample preparation. The formation 
energy, E/, is thus defined as 

Ef = £rot — EznO — ni[ii, (1) 

i 

where E^ ot is the total energy of the supercell with 
impurity, -EznO is the energy of pure ZnO host with 
the equivalent number of stoichiometric ZnO units, 



[ii is the chemical potential of a species and is 
the change of number of the corresponding species 
in the supercell. In thermodynamic equilibrium 
with bulk ZnO, the condition 

MZnO = A f o + MZn- (2) 

must always be satisfied. This leaves us with just 
one unknown, which we choose to be ^o- How- 
ever, an interval for the values of [io correspond- 
ing to physically meaningful growth conditions can 
be provided: The upper bound of the chemical po- 
tential, corresponding to O-rich growth conditions, 
is taken from the oxygen molecule, /io — \Eoi- 
Oxygen-poor conditions (or, equivalently, Zn-rich 
conditions) correspond to the equilibrium with 
metallic bulk Zn, /iz n = E^ lk , and hence, from 
Eq. ©, fxo = /IZnO - E^ k . 

A. Czn and Co impurities 

The formation energies of carbon at the Zn- 
site (Czn) and at the O-site (Co) obtained from 
calculations with a 3x3x3 supercell of ZnO are 
summarized in Table U for two different environ- 
ments. The limiting chemical potentials of the C- 
rich and for C-poor cases are taken as /ic = Eq 
and /ic = Eco — \ Eo 2 > respectively, as described 
in Ref . [2(3| • We note that the deviation of the val- 
ues as compared to Ref. [2(| is due to the different 
exchange-correlation functionals used in both cal- 
culations. The values obtained by Tan et at. (37j 
should not be directly compared to ours because 
they have used the total energies of single atoms as 
reference chemical potentials. In Fig. [2j the forma- 
tion energies which we obtain for Co and Cz n im- 
purities arc plotted as a function of oxygen chemi- 
cal potential. Our formation energies suggest that 
Co and Cz n impurities are energetically favorable 
only under carbon-rich conditions (negative values 
in the formation energy). Also, the formation en- 
ergy of Czn is smaller than the formation energy 
of Co in all chemical environments. 

We find that the formation of Co is energeti- 
cally favorable only in a narrow range of no shown 
as the thick (blue) vertical line. This tiny region 
corresponds to an equilibrium with metallic zinc, 
i.e., extremely O-poor conditions. Hence stabiliz- 
ing the Co defect in ZnO is difficult under com- 
mon experimental growth conditions of high oxy- 
gen pressure where /io approaches the value corre- 
sponding to molecular oxygen. 

A single Co in a ZnO supercell leads to a spin- 
polarized solution with an integral magnetic mo- 
ment of 2hb/C as discussed in the literature pre- 
viously [12, EH, HH . The spin polarized solution is 
lower in energy than the non-spin polarized one by 
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Table I: Formation energy of Co and Czn calculated from Eq. (JXJ) . The data show that the Cz n impurity is 
thermodynamically more stable than the Co impurity. 
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Figure 2: (Color online) The formation energy of Co 
and Czn defects as a function of no . The thick vertical 
line (blue colored line) shows the narrow range where 
the Co impurity can compete with Czn- 



a difference of 0.215 eV. The corresponding DOS of 
Co as single impurity in supercell S222 is shown 
in Fig. [3] (d), (e) and (f), employing the GGA, 
GGA+Ud and HSE06 functionals. For compar- 
ison, the DOS of pure ZnO obtained from each 
method is shown in Fig. [3] (a), (b) and (c). In 
GGA+Ud, the value Ud = 7.5 eV is applied to 
the (fully occupied) Zn ri-orbitals. We find a mag- 
netic solution for Co in all the methods. Note 
that both in GGA+U d and in HSE06, the impu- 
rity band in the minority spin is clearly separated 
from the valence band. The real-space plot of the 
effective spin density (from GGA) obtained by sub- 
tracting the charge densities of majority and mi- 
nority spins (Ap = p^ — Pi) is shown in Fig. 01 One 
finds that the effective spin density is centered on 
the C atom and extends to the nearest neighbor 
Zn atoms and the next-nearest neighbor O atoms 
with a total magnetic moment of 2/xs/C. We note 
that the HSE06 functional gives a large exchange 
splitting of the impurity states and larger band gap 
for ZnO than GGA, and the impurity states intro- 
duced by Co which mainly consist of carbon states 
are deep in the gap when measured from the va- 
lence band maximum (VBM). Therefore, it is not 
possible to use Co as an acceptor that could in- 
troduce hole states for thermally activated p-typc 
conductivity in ZnO. On the other hand, the more 
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Figure 3: (Color online) Density of states (DOS) of 
pure ZnO ((a), (b), (c)) and Co impurity in ZnO ((d), 
(e), (f)). The DOS in different rows are obtained from 
GGA, GGA+Ud on Zn and from HSE06 functional. 
The spin-up DOS and the spin-down DOS are plotted 
as positive and negative y-axis. The DOS of Co im- 
purity in ZnO in S222 geometry consists of the atomic 
contribution of nearest neighbor Zn, O to C in units 
of states/atom, while the total DOS is scaled to the 
number of formula units of ZnO. 



stable Czn impurity acts as a donor and does not 
lead to any spin polarization, in agreement with 
Ref. [H. 

Finally, we discuss possible compensation effects 
that could alter the magnetic properties of Co- 
Using DFT studies, Li et al. have shown that 
the presence of oxygen vacancy ( Vo ) quenches the 
spin-polarization from Co in ZnO [40j. Compar- 
ing the formation energy of neutral defects [4l| . 
we conclude that both, Vo and Vz n , have positive 
formation energy, with the formation energy of Vo 
being lower than the formation energy of Vz n in 
O-poor conditions (favorable conditions for Co as 
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Figure 4: (Color online) Effective spin density (Ap = 
Pt ~ Pi) plotted along the plane perpendicular to c- 
axis and passing through Co- The light (magenta), 
dark (turquoise) and black colored balls are used to 
represent Zn, O and C atoms [42]. The isolines (and 
also the color gradient) are plotted from 0.003 eA 
to 0.028 eA~ 3 at intervals of 0.005 eA -3 . One can 
clearly observe that the spin-polarization arising from 
Co extends to relatively large spatial region. 
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Figure 5: (Color online) DOS of Co in supercells S122 
(a), S221 (b) and S222 (c), respectively. The DOS 
consists of the atomic contribution of nearest neighbor 
Zn and O to C (in units of states/atom), while the 
total DOS is scaled to the number of formula units of 
ZnO. 



seen above). Hence, the role of defects like Vz n 
and Vo can be safely ruled out in contributing to 
magnetism in the material. However, the magnetic 
moment of Co can be quenched by H impurities in 
ZnO that act as donor. We have tested co-doping 
of H and C in S622 and find that the electron from 
H compensates one Bohr magneton of magnetic 
moment. Upon adding two H in the supercell, the 
magnetic moment is completely quenched to zero. 
This observation is independent of the distance of 
the hydrogen from the carbon site in the super- 
cell. In experiments the hydrogen concentration 
during ZnO synthesis is difficult to control. It has 
been found that the energy barrier for the- escape 
of H through the ZnO surface can be as large as 
0.58 eV [43|. Thus even if one succeeds with Co 
doping, one may still fail to obtain magnetic mo- 
ments due to the ubiquitous presence of hydrogen. 
Similar conclusions have also been derived from 
formation energy studies [H| . 



B. Fixed spin moment calculations using 
supercell approach 

Although the thermodynamic stability consid- 
erations of the previous section point to the dif- 
ficulty in incorporating larger amounts of carbon 
on the oxygen site in ZnO, we continue to discuss 
the consequences of Co impurities for the magnetic 
properties, assuming that an appreciable concen- 
tration can be built up by non-equilibrium prepa- 
ration techniques. In order to investigate impurity 
interactions at various concentrations, we perform 
DFT calculations with both VASP and FPLO us- 



ing a single Co impurity in the different supercells 
shown in Fig. [T] (a) to (c) . These cells correspond 
to carbon concentrations of 12.5% (for both S122 
and S221) and 6.25% (S222). With decreasing the 
C concentration from to 12.5% and 6.25%, the in- 
teraction between impurities becomes weaker and 
the impurity band becomes narrower. Even though 
the supercells S221 and S122 have the same doping 
concentration of 12.5%, they do not lead to equiva- 
lent electronic structure because of the anisotropy 
of the wurtzite crystal. From Table [U we note 
that the magnetic moment of S122 is a fractional 
value in units of hb, an indication of strong hy- 
bridization between the C 2p orbitals and the va- 
lence band, leading to metallic behavior as shown 
in Fig. [5ta). However, the magnetic moment of 
S221 and S222 is integer (2 /^b), i.e., the system 
is a magnetic half metal (see Fig. [5] (b) and (c)). 
We carried out fixed-spin-moment (FSM) calcula- 
tions to test the stability and characteristics of the 
spin-polarized solutions itself. The plot in Fig. [6] 
shows the magnetic polarization energy of the su- 
percells shown in Fig. [1] (a) to (c). Here the mag- 
netic polarization energy is defined as the energy of 
the system at a given fixed magnetic moment (m) 
relative to the energy of non-magnetic solutions 
(m = 0). When a spin-polarized solution is fa- 
vored, then the magnetic polarization energy takes 
negative values. For metallic systems (supercell 
S122), the minimum magnetic polarization ener- 
gies are found at non-integer values of m as shown 
in Table HH and the magnetic polarization energies 
show a smooth parabolic change around the min- 
ima. For supercells S221 and S222 the minimum 
magnetic polarization energy is obtained at an in- 



6 



Table II: Nearest neighbor (nn) and next-nearest neighbor (mm) distance of the periodic image of the C impurity 
in A, magnetic moments as obtained from VASP (FPLO), and formation energy (£7) for O-poor and C-rich 
conditions. The total magnetic moment (M-rotai) is shown per supercell, while the magnetic moments of Zn 
(Mzn) and O (Mo) are shown per atom, averaged over the nearest neighbor Zn and O, respectively. 
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Figure 6: (Color online) Magnetic polarization energy 
(E ap — E nsp ) (sp = spin polarized and nsp = non- 
spin polarized) as a function of the magnetic moment 
for different supercells from fixed-spin-moment (FSM) 
calculations. The minima in the magnetic polariza- 
tion energy indicate the most stable magnetic moment 
in the supercell. The magnetic polarization energies 
obtained from VASP and that of FPLO are in close 
agreement with each other. 



teger value m = 2 /xb and the energy has a cusp at 
the minima which correspond to localized impurity 
band in the gap. Thus, from the FSM calculations 
we verify the stability and characteristics of spin- 
polarized solutions in the supercells and one can 
estimate an upper bound of C concentration for 
the material to be a semiconductor rather than a 
metal. A higher C concentration leads to a metal- 
lic solution, i.e. the impurity band hybridizes with 
the valence band of ZnO. According to the GGA 
results, the distance between C impurities should 
be larger than twice the lattice constant, meaning 
that the C concentration should be lower than 6%, 
to avoid such a metallic state. The precise esti- 
mate of this concentration might depend on the 
exchange-correlation functional used in the calcu- 
lations. As we discussed above, hybrid-functionals 
predict the impurity bands to lie further up in 
the gap. Consequently, ZnO may tolerate an even 
larger C concentration before becoming metallic 
than the values we have estimated from GGA. 



C. Exchange interaction in C doped ZnO 

Next, we study the magnetic interaction be- 
tween the carbon impurities as a function of their 
separation. We employ the relatively large super- 
cells S622 and S226 (Fig.ffld and e) to include long- 
range interactions. The concentration of carbon in 
these supercells is 4.17% which is within the range 
of experimental impurity concentrations [2l| . As 
long as the carbon impurities are well-separated 
(maximum separation of 16 A in c-direction and 
11 A in a6-direction for S226 and S622, respec- 
tively), each impurity retains its magnetic moment 
of 2\ib independent of the distance. Therefore, it 
is reasonable to map the interaction between two 
Co impurities onto an effective Ising Hamiltonian 

H(r i0 ) = -- ^2 Jio{no)(Ti(To . (3) 

(i) 

The index i runs over the nearest neighbor sites 
appropriately taking into account the periodic 
boundary conditions used in the DFT calculations. 
Without loss of generality we can choose gq = ±1 
and tTj = 1 to calculate the individual interaction. 
For a given geometry — \n — ro|, the J^o can 
now be accessed from the total energy difference 
between the ferromagnetic and antifcrromagnctic 
alignment of the moments of Co impurities. 

These J;o are plotted in Fig. [7J As a guide to the 
eye, we include a shaded area bounded by two ex- 
ponentially decaying functions within which all the 
interactions arc contained. We note that the Jio 
are positive in the whole range, i.e. the spin inter- 
action between Co is ferromagnetic in the host ma- 
terial. Averagely, the decrease with increasing 
separation r^o between the impurities while there is 
a large spread in the data for small distances. We 
categorize the geometric configuration of the impu- 
rities using the crystallographic direction and se- 
lect those configurations that have the major com- 
ponent of their distance vectors in the hexagonal 
plane (afc-direction). The selected configurations 
(Al, A3, A4, A6 and A12) are shown as filled sym- 
bols in Fig. [Jj We observe a strongly anisotropic 
behavior of J;o; the filled symbols are located in 
the upper part of the shaded area, whereas the 
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Figure 7: (Color online) Exchange interaction ( Jio) ex- 
tracted from the total energy using Eq. ((3]) plotted 
against the relaxed Co-Co distance (r;n). CI - C13 
and Al - A13 represent the different positions of Co 
in supercell S226 and S622 as in Fig. Q] (d) and (e), 
with the other Co positioned at CO and AO, respec- 
tively. The range of interaction values is bounded by 
two exponential curves. The Jin for the Co-Co ori- 
ented along the afo-hexagonal plane, which are shown 
as filled symbols, are larger than those along the c-axis. 
From the minority-spin band structure of a Co impu- 
rity in ZnO (in CCA) shown in the inset one finds that 
the Co band crosses the Fermi energy in the hexagonal 
plane (A-L, T-M, A-H and F-K directions) imparting 
metallic character and thus the magnetic interactions 
along the hexagonal plane are stronger. 



open symbols are in the lower part. Test calcula- 
tions performed for some configurations with the 
GGA+Ud method have confirmed the ferromag- 
netic nature of the interactions and indicate that 
the results are independent of the exact energetic 
position of the minority-spin C 2p-states within the 
host band gap. 

Since the carbon impurity levels lie deep in- 
side the band gap (see the DOS in Fig. [3] (d)), 
the magnetic interactions cannot be mediated by 
carriers in the valence or conduction bands of 
the host material (as typical for small-gap DMS 
materials). Instead, ferromagnetism arises from 
partly filled impurity bands, as in Zener's double- 
exchange model [lj, HH, Hg. According to this 
model, ferromagnetic ordering becomes stabilized 
because it allows for a lowering of the kinetic en- 
ergy of the electrons in the impurity band: Due to 
the overlap of impurity orbitals of the same spin, 
the dispersion of the impurity bands is increased, 
and the kinetic energy of electrons occupying the 
lower part of these bands is decreased. Obviously, 
this energy gain is largest for interactions along 
those crystallographic directions where the impu- 
rity band crosses the Fermi level. 

The anisotropic nature of the magnetic inter- 



action can therefore be explained from the band 
structure of Co in ZnO. In the inset of Fig. the 
band structure of the minority spin states of Co 
in S333 is shown. Since the spin-majority impu- 
rity bands are completely filled and thus cannot 
give rise to any energy gain, only the minority-spin 
bands can be responsible for the magnetic interac- 
tions. We find that the carbon bands are cross- 
ing the Fermi level mostly along the hexagonal ab- 
plane, i.e., along the A-L, T-M, A-H and T-K 
directions, while along the c-axis (M-L, T-A and 
K-H) the Fermi energy lies in the gap between the 
occupied and the unoccupied bands. Therefore, 
the ferromagnetic interaction is mediated more ef- 
fectively within the hexagonal a6-planc, giving rise 
to anisotropic ferromagnetism. 

We note that the ferromagnetic interaction is 
rather short ranged (below 10 A). Ferromagnetism 
at finite temperatures can thus only be expected 
at a carbon concentration larger than 2% (refer to 
the scale of the C concentration in the upper axis 
of Fig. [7]) . Together with the estimate given in Sec- 
tion IIIIB1 which gives the upper bound for C con- 
centration to avoid metallic state, we conclude that 
a Co concentration between 2% and 6% should be 
optimal to obtain ferromagnetism due to homoge- 
neously distributed impurities. Moreover, we note 
that C impurities have a tendency to cluster which 
poses an upper limit to the useful concentration. 
As pointed out in previous work [47] , the C im- 
purities tend to form energetically more stable C2 
molecules. In this case, the ppn* orbital of the 
C2 molecule resonates with the conduction band 
of ZnO which may lead to ferromagnetism medi- 
ated by the host material in n-type ZnO. Thus, 
this mechanism differs from the ferromagnetic in- 
teraction between Co impurities described in the 
present work which originates from conductivity 
due to impurity bands, and is hence independent 
of host carriers. 



IV. CONCLUSIONS 

We have studied the two possible substitutional 
carbon impurities in ZnO (Co and Cz n ) and find 
that Czn is energetically more favorable than Co- 
However, Co can be stabilized under specific (O- 
poor and C-rich) growth conditions. This type of 
environment is uncommon in usual ZnO growth, 
but could be generated by oxidizing metallic zinc 
in the presence of an atomic carbon source. Given 
that a material with a Co concentration of 2-6% 
can be realized, this would be an interesting pro- 
totypical system for d° magnetism: Our DFT cal- 
culations indicate that Co in this range of con- 
centrations is associated with a localized magnetic 
moment of 2 (J-b/C. Partially filled impurity bands 
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in the minority spin channel mediate ferromagnetic 
interaction between the Co impurities. These in- 
teractions are short-ranged and anisotropic, being 
stronger within the a6-plane of the wurtzite ZnO 
crystal than along the c-axis. Based on our calcula- 
tions, we predict that such layered fcrromagnetism 
could be used as an experimental hallmark to dis- 
tinguish carbon-induced magnetism from possible 
other forms of defect-induced fcrromagnetism in 
ZnO. 
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